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Abstract 

The time-evolution operator for an explicitly time-dependent Hamiltonian is ex- 
pressed as the product of a sequence of unitary operators. These are obtained by 
successive time-dependent unitary transformations of the Hilbert space followed by 
the adiabatic approximation at each step. The resulting adiabatic product expansion 
yields a generalization of the quantum adiabatic approximation. Furthermore, it leads 
to an infinite class of exactly solvable models. 

PACS number: 03.65.Bz 

Consider the dynamics of a quantum mechanical system whose Hamiltonian H = H{r) 
is explicitly time-dependent. The evolution of a state vector \ip(r)) is governed by the 
Schrodinger equation: 

ih\i){r)) = H(t) |V(t)) , |V(0)) = |-0o> , (1) 

where a dot means a time-derivative. Alternatively, one has \ip(r)) = U(t)\i/;o), where 

U(t) = T e-^^ H{t)dt (2) 

is the time-evolution operator. Here T denotes the time-ordering operator. The purpose 
of this article is to express U(r) as the product of a sequence of unitary operators each of 
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which is the adiabatically approximate time-evolution operator in some adiabatically moving 
frame. 

In order to construct these unitary operators, let us separate the adiabatic part U^°'(r) 
of the exact time-evolution operator: 

U(t) = U w (t) V {1 \t) . (3) 

U^(t) is defined by 

U (0) {r) := JV Q " (r) |n;r)(n;0| , (4) 

n 

a n (r) := 5 n (r)+7„(r), (5) 

1 f T f T d 

$n(r) := -- / E n (t)dt, 7 n (r) := i \ (n;t\ — \n;t)dt, (6) 
h Jo Jo at 

where \n;t) are instantaneous eigenvectors of the Hamiltonian H(t) with eigenvalue E n (t), 
i.e., 

H(t)\n;t)=E n (t)\n;t) . (7) 

Throughout this article, it is assumed that the energy spectrum is discrete, all the eigenvalues 
of the Hamiltonian are non-degenerate and there is no level-crossing. 

The adiabatic time-evolution operator U^(t) as defined by (f|) includes the effects of 
Berry's geometric phase fTJ, through the phase angles 7 n (t)- The operator V^(r) includes 
the non-adiabatic effects. As U^(t) is unitary, Eq. (§) may be viewed as the definition of 

Next consider an arbitrary time-dependent unitary transformation of the state vectors: 
\iP{t)) — > \iP'(t)) := U(t)\iP(t)) . This may be viewed as transforming to a "moving frame" 
of reference. In the "moving frame" the state vectors also satisfy a Schrodinger equation 
with a new Hamiltonian T~C(t). This is related to H{r) according to 

H{t) = U{t) H{t) U\t) - %h U{t) U\t) . (8) 

Now let us set U(t) = U^(r). Then one can easily check that the state vectors |^'(t)) 
in the "adiabatic moving frame" evolve according to W{t)) = V^(r)\^'(0)) = V^(r)\ip ). 
Thus, in view of Eq. (j8|), we can express V^'(r) directly in terms of the corresponding 
transformed Hamiltonian which we shall denote by H^(r), namely: 

l/W(r)=Te-^ H(1,(1)di . (9) 
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Since V^(t) is also a time-evolution operator, the same procedure may be repeated for 
V^\t). In other words, we can use (r) in place of the original Hamiltonian H(t) 
to define the adiabatic part U^(r) of the time-evolution operator V^{r) and therefore 
write VW(t) = U^ 1 \t)V^(t). Clearly this can be continued indefinitely. The result is an 
adiabatic product expansion of the time-evolution operator: 

oo 

U(t) = UU®(t). (10) 

i=0 

This expression may be viewed as a generalization of the quantum adiabatic approxima- 
tion This generalized quantum adiabatic approximation is carried out by calculating the 
first N terms in the product in ( [TDD an d approximating U(t) by 

TV 
i=0 

Incidentally the transformed Hamiltonian H^\t) may be easily obtained in the eigenbasis 
of the initial Hamiltonian. This is done by substituting U^°"(t) for U{r) in ([8]). The result 
is 

#W(r) = -%% e- i[Qm(r) " a " (r)1 A mn (r) |m;0)(n;0| , 

where A mn (r) := (m; r|-r-|n; r), and the second equality is obtained by differentiating Eq. (]?]) 
with respect to time and using orthonormality of the energy eigenvectors. The appearance of 
the time-derivative of the Hamiltonian in the expression for #«(t) and the fact that #«(r) 
is off-diagonal in the eigenbasis of H(0) are reminiscent of the meaning of the adiabatic 
approximation in which (r) is neglected. 

Eq. (|T2| ) can also be used to yield the transformed Hamiltonian H (r) corresponding 
to the (i + l)-th term in (TTOj) . This is done by introducing the symbol nS % > which labels 
the eigenvectors and eigenvalues of H^(r), and replacing H^(r), n, m, H(t) in (0) by 
(r), Tny\ and H^'(t), respectively. Clearly, H®(t) and consequently 

U®(t) :=Y,e ia ^ {T) |n (i) ;r)(n w ;0| , (13) 



3 



involve i-th time derivatives of the original Hamiltonian. 

The generalized quantum adiabatic approximation (|TT| ) of order N is a reliable ap- 
proximation for the time-evolution operator, if one can neg lect H( n+1 \t). It is exact, if 
H^ n+1 ^(t) = 0. In fact, the latter equation may be used as a defining condition for gener- 
ating a class of exactly solvable examples for which U(r) = U^(r) ■ ■ • U^ N \r) is the exact 
time-evolution operator. 

In the remainder of this article, I shall demonstrate the application of the adiabatic prod- 
uct expansion in the analysis of the dynamics of a magnetic dipole in a changing magnetic 
field. 

The Hamiltonian of this system is given by 

H[R] = H{r,9,<p) = bR{r,6,(f) ■ J = br (sin 9 cos <p.J l + sin ^ sin J 2 + cosflJ 3 ) , (14) 

where b is the Larmor frequency, (r, 9, (p) are spherical coordinates and J is the angular 
momentum operator with components J M , /i = 1,2,3. The condition r / guarantees 
the lack of level crossing. Therefore the appropriate parameter space M of this system is 
IR 3 — {0}. The time-dependence of the Hamiltonian corresponds to the curve C : [0, T] — > M 
traced by the tip of the magnetic field in time. This curve is parametrically given by 
(r(t),9(t),<p(t)). Assuming that for all t £ [0, T], C(t) ^ (r,9 = ir,(p), one can choose a 
complete orthonormal set of single-valued eigenvectors \n; R) of H[R\. These are given by 
I: 

|n; (r, 9, <p)) = |n; (r , 9,<p)) = W(9, V )\n; (r , 0, 0)) , 8 £ [0, tt), <p £ [0, 2vr) , (15) 
where 

W(9, if) := e-^ j3 e-f j2 e^ j3 , (r , 9 , ^ ) := (r(0), 9(0), <p(0)). (16) 
The corresponding eigenvalues are: 

E n (r, 9, <f) = E n {r, 0,0) = h nbr , with n = 0, ±^ , ±1 , ±^ , • • • . (17) 

Note that in (|T3j), \n; (r, 0, 0)) are the eigenvectors of H(r, 9 = 0, (p = 0) = brJ 3 , i.e., 

J 3 \n; (r,0,0)) = hn\n; (r,0,0)) . 

Assuming that the eigenvectors \n; (r,9,tp)) of the Hamiltonian are also eigenvectors of 
| J| 2 , i.e., they have definite angular momentum: 

|J| 2 \n;{r,9,ip)) = ](] + !) \n; {r,9,ip)) , n = -j, -j + 1, • • • ,j , (18) 
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one can calculate 



A mn (t) 

A(mn) 
j^rnn) 



4 (ran) 



5 n (t) 



A (mn) ■ (t) + A (ran) ^ + A ( mn) ^ 



9 

(m; (r,0,<p)| — |n; (r,0,<p)) 







<9 1 . 

(m;(r,0,</?)| — |n;(r,0,</?)) = - (e t¥, C m <5 mn -i - e ^C^ra-in) 

(m; (r,9,(p)\— \n; (r,6,<p)) , 
m(l - cos Q)5mn + g sin61 (^Cm^mn-i + e~ %<p C n 6 m - lTI 

w(9 , (po ) \m j + + n*(t) j-} w*(e ,<p ) , 



-ti6 / r(0 dt' 

JO 



7n (t) = -n7(t) 



where C m := w (j — m)(j — m + 1), W is defined in Eq. (0), 



(19) 
(20) 

(21) 

(22) 
(23) 
(24) 



;= e -*K*)+v(*)l [ S m9(t)tp(t)+i9(t)), a(t) := 5(t) + 7 (t) , (25) 
5(t) : = -b [ l r{f) dt' , 7(t) := - cos 6(t')]p(t')dt' , (26) 

and extensive use is made of the properties of J^ 1 and := J 1 ±iJ 2 , particularly: 

J± |m; (r, 0, 0)) = AlC ±m \m ± 1; (r, 0, 0)) . 

), and carrying out the necessary computations, one 



Substituting (U) and (H) in 
obtains: 



/ \ I r 9 ^0 / \ *? ^0 / / \ \ i 1 

Ct;(i) < [cos — coscr(t) — sin — cos(2</?o + cr (t))J J + 



[— cos 2 — sin a{t) — sin 2 ^ sin(2<^ + cr(t))) J 2 + [— sin 9 cos cr(£)] J A 



0o 



(27) 



: 6 r (1) (t) [sin (1) (t) cos <^ {1) (t) J 1 + sin (1) (t) sin <^ (1) (t) J 2 + cos (1) (t) J 3 1 , (2 



where u;(£) := |fi(t)| — y9 2 + sin 2 9 (p 2 is the angular speed of the tip of the magnetic field, 
a(t) is the phase of fl(t), i.e., 



cr(i) := —a — ip + £ mod 27T . 



cos£ := 



sin 9 (p 



■ t 9 
sm£ := - 



and (t) , (t) , ipW (t) are defined by equating the right hand sides of ( p7[) and fl28|), i.e., 

u(t)A(t) 



r (1) (t) := 



cos 



tan 



:= yl + sin^osin #0 sin[<£> + 2cx(i)] 



sin 9 cos cr(t) 

AW , 
sin a(t) + tan 2 (f ) sin[2y2 + a(t)] 



(29) 
(30) 

(31) 



-coscr(t) +tan 2 (| i )cos[2v9 + a(t)]_ 

This notation allows one to directly read off the eigenvectors In' 1 '; t) and eigenvalues E n w (t) 
of H^'(t) by replacing n by and (r,9,p) by (r^, y^ 1 -*) in Eqs. ( p~5| ) and (|T7|). Par- 
ticularly interesting is the fact that r^'(t) > 0, unless the angular speed of the tip of the 
magnetic field vanishes. u(t) = = r^'(t) define the times t for which H^(t) = 0. These 
"degenerate" situations need further treatment. Here I suffice to remind the reader that for 
time periods during which H^'(t) = 0, the adiabatic approximation is exact. 

Eqs. (H)-([DD together with ©, (0), and (g) can be used to obtain U^ i+1) (r) for all 
% = 0, 1, • • -. This is done by substituting r«, 0«, for r, 9, cp in Eqs. ©-(0), (HI), 
(0)) (S) an d using (|i3|). The result is: 

[/(*+%) = w(0«(T)^ (l ' ) (^))E eia " w(T V ) (ro,O,O))(n,(r Q ,O, 0)|W (0«(O), ¥>«(())) , 



where a w (t) := 5 {i) {t) +7 (i) (t) and 



e*^- 73 iy (0«(O),p (i) (O) 



5 (l) (t) 
a„(i) (£) 



-6 / r®{t')dt', 7 w (t) := 



[1 - cos 9^ (t')} p^(t') dt' . 



na®(t), 5 n(i) (t) = n5®(t), 7nW (t) = n 7 (i) (0- 



For yj = Eqs. (p9|)-(|3T1) simplify considerably. In this case one has: 



r (1) (t) 



#W (t) = cos' 1 [- sin 9 cos a(t)} , (t) = - tan" 1 



tancr(t) 
cos 9n 



(32) 

(33) 
(34) 

(35) 



Substituting these relations in Eq. ( |33"D with i = 1, and performing the necessary algebra, 
one finds: 



*W(t) 



v(t')dt' =: -£(t) , 7 (1) (t) 



tan 



-i 



x(t) + y(t) 

l-X(t)F(t) 



(36) 
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, w . cos^nftancrft) — tancr(0)l . „ r . . , . . 

Xt := ° L — 1 \, u 77^, y * := tan 5 sin at -sinaO , 37 

cos^ Vq + tan a{t) tan cr(0) 

where £(t) is the length of the projection of the portion of the curve C which lies between 
C(0) and C(t) onto the unit sphere centered at r = 0. The appearance of such a geometric 
quantity in this calculation is quite remarkable. 

Let us next apply the adiabatic product expansion to identify a class of exactly solvable 
cases. This is done by enforcing H^ n+1 \t) = or alternatively r^ N+1 \r) = 0. The first 
nontrivial case is N = 1, i.e., r^(r) = 0. This condition may be easily fulfilled by requiring 
a = which immediately leads to: 



r(0 = I 



cos 8 ip — 



dt v sin 9 ip 1 



b 



d r 8' 

n dip v sin t 
COS 6 ; 



1 + 



(38) 



v sm o ip' 

where 9' = d9/d(p. Thus, for every 9 = 6(<p) which renders the right hand side of (38) positive 
for all t, \ip(t)) = U^°> (t)!!^ 1 ' (t)\^o) is the exact solution of the Schrodinger equation ([!]). 
Clearly, one can generate further exactly solvable examples corresponding to N — 2, 3, • • •, 
in this manner. 

The adiabatic product expansion suggested in this article may be viewed as an iterative 
procedure to integrate the Schrodinger equation. It involves only algebraic manipulations 
and integration of scalar functions. It also provides a generalization of the quantum adiabatic 
approximation which involves an integer N as its order. The approximation yields the exact 
result if H^ N+1 \r) = 0. In the general case where H^(r) 0, for all i, one may improve 
the accuracy of the approximation indefinitely by increasing its order. 

As shown for the magnetic dipole system, one may be able to obtain generic expressions 
for the terms in the adiabatic product expansion. These can be used to generate exactly 
solvable examples. In particular, if the dynamics of the system is given by a dynamical 
group G, i.e., if the Hilbert space provides a unitary irreducible representation p of a Lie 
group G and H(t) = p[h(t)}^ H d (t)p[h(t)}, with h(t) G G and H d (t) diagonal, belongs to 
the representation of the Lie algebra of G and U(t) = p[g(t)] for some g(t) G G [£|, then 
one might try to employ a similar approach as the one used in the analysis of the dipole 
system (G = SU(2)). In particular this approach may be readily applied for the generalized 
harmonic oscillator system where G = SU(1,1), ||. An interesting direction for further 
study of the adiabatic product expansion is to seek its meaning and implications in the path 
integral approach to quantum mechanics and ultimately quantum field theory. 



7 



Acknowledgements 

I would like to thank Bahman Darian and Ertugrul Demircan for reading the first draft of 
this article. I would also like to acknowledge the financial support of the Killam Foundation 
of Canada. 

References 

[1] M. V. Berry, Proc. Roy. Soc. London A392, 45 (1984). 

[2] A. Messiah, Quantum Mechanics, Vol. 2, North-Holland, Amsterdam, 1962. 

[3] A. Bohm, Quantum Mechanics: Foundations and Applications, third ed., chapter 22, 
Springer- Verlag, New York ,1993. 

[4] A. Mostafazadeh, J. Math. Phys. 37, 1218 (1996), and references therein. 

[5] R. Jackiw, Int. J. Mod. Phys. A2, 285 (1988), and references therein. 



8 



